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COMPUTER  GENERATED  HOLOGRAPHIC  OPTICS 


STATEMENT  OF  THE  PROBLEM  STUDIED 

The  goal  of  this  research  was  to  theoretically  and  experimentally 
Investigate  the  use  of  computer-generated  holograms  In  the  design  and  fabri¬ 
cation  of  holographic  optical  elements  (holographic  lenses  or  diffraction 
optics).  Conventional  design  of  holographic  lenses  had  previously  been 
limited  to  spherical  recording  beams,  and  more  recently  to  recording  beams 
formed  by  conventional  optical  systems.  The  utilization  of  computer-generated 
holograms  In  the  recording  beams  allows  a  completely  arbitrary  specification 
of  the  recording  wavefronts  in  the  hologram  plane,  making  possible  a  hologram 
lens  of  higher  performance  (lower  aberrations)  than  previously  possible. 

We  refer  to  an  interferometrlcally  recorded  hologram  lens  having  a  computer¬ 
generated  hologram  In  a  recording  beam  as  a  computer-originated  holographic 
lens.  We  studied  the  basic  limitations  of  computer-generated  holograms,  and 
integrated  the  analysis  of  computer-generated  holograms  and  computer-originated 
holographic  optical  elements  Into  ERIM's  holographic  optical  analysis  and 
design  (HOAD)  ray-tracing  computer  program.  In  addition,  we  designed, 
fabricated,  and  tested  an  optical  system  using  a  computer-generated  hologram 
and  a  computer-originated  holographic  lens. 

SUMMARY  OF  RESULTS 

Several  different  types  of  computer-generated  holograms  were  studied  and 
evaluated  for  t hi s  application.  Important  considerations  are  the  diffraction 
efficiency,  the  type  of  optical  recording  device  required  (binary  vs  grey- 
level),  the  spatial  efficiency  with  which  the  recorder  Is  utilized,  the 
type  of  recording  material  required  (absorbing  vs  phase  shifting),  the 
recording  material  properties,  and  the  Intrinsic  errors  (e.g.  spurious 
diffracted  terms)  arising  from  the  computer-generated  hologram.  For  achieving 
high  performance  holographic  lenses,  the  most  Important  consideration  of 
those  listed  above  Is  the  minimization  of  Intrinsic  errors  and  spurious 


terms.  Given  the  availability  of  a  grey-level  recording  device,  we 
concluded  that  the  simple  carrier  method  is  the  computer-generated  hologram 
best  suited  for  this  application.  It  consists  of  a  carrier-frequency  en¬ 
coding  of  the  phase  plus  a  bias,  resulting  in  a  real,  nonnegative  transmittance 
of  the  form:  bias  +  amplitude  •  cos  [wx  +  phase  (x,  y)]. 

Computer-generated  holograms  are  not  very  useful  as  optical  elements 
by  themselves  since  they  have  low  diffraction  efficiency,  low  numerical 
aperture  (i.e.  low  space-bandwidth  product),  and  intrinsic  spurious  orders 
of  diffraction.  However,  all  of  these  deficiencies  can  be  avoided  by  inter- 
ferometrically  recording  a  holographic  optical  element  on  a  volume-phase 
material  of  high  diffraction  efficiency,  such  as  dichromated  gelatin,  using 
a  computer-generated  hologram  in  one  of  the  recording  beams.  Referred  to  as 
a  computer-originated  holographic  optical  element  (COHOE),  it  can  have  a  high 
numerical  aperture  if  one  or  both  of  its  recording  beams  is  produced  by  a 
conventional  optical  system  in  conjunction  with  the  computer-generated 
hologram.  This  is  described  further  in  Appendix  A.  Appendix  A  also  describes 
the  Implementation  of  an  analysis  capability  for  computer-generated  holograms 
and  COHOEs  within  the  HOAD  program.  The  design,  analysis,  fabrication  and 
testing  of  an  optical  system  employing  a  computer-originated  holographic 
lens,  using  a  computer-generated  hologram  in  one  of  Its  recording  beams* 
is  also  described.  By  this  method  we  arrived  at  a  holographic  lens  having 
far  lower  aberrations  than  a  holographic  lens  made  with  conventional  spheri¬ 
cal  recording  beams.  When  used  as  a  Fouier  transform  element,  this  aspheric 
holographic  lens  achieves  close  to  an  order  of  magnitude  improvement  in  two- 
dimensional  space-bandwidth  product  capacity  as  compared  with  a  conventionally 
recorded  holographic  Fourier  transform  element.  This  new  technology,  which  we 
have  demonstrated  to  be  feasible,  Is  expected  to  generally  provide  better 
performance  with  fewer  optical  elements  than  what  was  previously  obtainable. 

In  the  course  of  examining  the  Imaging  of  computer-generated  holograms 
from  one  plane  to  another,  we  studied  the  Imaging  of  wavefronts.  In  contrast 
to  imaging  an  Intensity  distribution,  the  Imaging  of  a  wavefront  by  an  optical 
system  Is  complicated  by  the  fact  that  the  phase  of  the  wavefront  In  addition 
to  Its  intensity  must  be  handled  properly.  Except  for  special  cases,  there  is 
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an  additional  quadratic  phase  factor  present.  Analysis  was  performed  to 
compute  the  wavefront  image  position,  magnification  and  additional  quadratic 
phase  factor  for  general  one-lens  and  two-lens  imaging  systems.  It  was  found 
that  by  using  two  lenses  of  fixed  focal  length  It  Is  possible  to  obtain  both 
a  given  magnification  and  a  given  quadratic  phase  factor,  within  certain 
limits.  The  range  of  obtainable  magnifications  and  quadratic  phase  factors, 
it  was  also  found,  can  be  expanded  by  using  three-lens  or  four-lens  Imaging 
systems. 

Application  of  computer-originated  holographic  optics  to  variable  phase 
compensators  was  also  analyzed.  It  was  determined  that  three  holographic 
phase  plates  can  be  used  to  produce  variable  (adaptive)  amounts  of  one¬ 
dimensional  cubic  and  quartic  phase  correction  by  the  translation  of  the 
second  and  third  phase  plates  with  respect  to  the  first  phase  plate.  Used 
in  conjunction  with  a  conventional  optical  processor  for  synthetic  aperture 
radar  imaging, these  phase  plates  would  allow  for  the  compensation  of  several 
wavelengths  of  phase  error  in  the  radar  signal  history.  This  task  is  current 
ly  very  difficult  to  perform  in  a  dynamic  manner  using  conventional  optics. 

Computer  experiments  and  analysis  were  also  performed  on  a  special 
holographic  diffuser  that  would  cause  a  minimum  of  speckle  in  a  reconstruct¬ 
ed  image.  The  holographic  diffuser  would  be  made  from  a  computer-generated 
hologram.  It  is  briefly  described  in  Appendix  B. 
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Abstract 


The  capability  of  analyzing  hologram  lenses  recorded  with  arbitrary  wavefronts  has  been 
added  to  a  holographic  raytracing  design  program.  The  recording  wavefronts  are  defined  by 
analytical  phase  functions,  for  example,  a  two-dimensional  polynominal  expansion.  The  co¬ 
efficients  of  the  functional  representations  of  the  recording  wavefronts  are  used  as  para¬ 
meters  to  optimize  the  performance  of  an  optical  system  containing  the  hologram  lens.  The 
optimum  recording  wavefronts  are  then  produced  with  the  help  of  computer-generated  holo¬ 
grams.  Several  useful  arbitrary  wavefront  phase  functions  are  discussed.  Design  predict¬ 
ions  and  experimental  results  will  be  shown  for  a  holographic  Fourier  transform  lens 
recorded  with  the  aid  of  a  computer-generated  hologram. 

Introduction 


The  past  several  years  has  seen  an  increased  use  of  refractive  optical  elements  having 
aspheric  surfaces.  The  generalization  of  lens  surface  profiles  has  produced  better  system 
performance  with  fewer  elements  in  many  applications.  It  is  reasonable  to  suppose,  there¬ 
fore,  that  holographic  optical  elements  (HOEs)  would  also  benefit  from  the  use  of  general¬ 
ized  recording  wavefronts  (Figure  1) .  The  design  of  diffractive  optical  systems  has  until 
recently  been  restricted  to  the  use  of  HOEs  recorded  with  plane  and/or  spherical  wavefronts 
(Figure  2) .  These  wavefronts  are  the  ones  most  readily  generated  in  the  laboratory  using 
standard  refractive  optical  components.  The  generation  of  truely  arbitrary  recording  wave- 
fronts  using  refractive  optics,  on  the  other  hand,  is  difficult  at  best.  A  more  feasible 
approach  is  the  use  of  a  computer-generated  hologram  (CGH)  to  optically  generate  the  de¬ 
sired  arbitrary  aspheric  recording  wavefront1.  The  increasing  availability  and  perform¬ 
ance  of  devices  for  recording  CGHs  makes  this  an  attractive  approach. 

This  paper  describes  the  implementation  of  an  aspheric  HOE  design  capability  within  an 
existing  holographic  raytrace  program.  Consideration  is  given  to  methods  for  defining  ar¬ 
bitrary  recording  wavefronts  on  curved  as  well  as  flat  substrates.  In  addition,  techniques 
for  defining  CGHs  to  generate  the  desired  arbitrary  recording  wavefronts  are  presented. 
Finally,  the  design  of  a  simple  aspheric  HOE  to  be  used  as  a  Fourier  transform  (F.T. )  lens 
is  described  and  experimentally  evaluated. 


Figure  1.  Aspheric  HOEs  are  recorded  with  an 
object  wavefront  (a)  derived  from  an  auxiliary 
optical  system  or  (b)  defined  analytically. 


Figure  2.  Conventional  (spher 
ical)  HOEs  are  recorded  with 
spherical  and/or  plane  wave- 
fronts  . 


Raytracing  through  HOEs 


A  short  discussion  of  the  raytrace  grating  equations  is  appropriate  at  this  point.  We 
assume  that  a  HOE  is  recorded  with  an  object  beam  (Obj)  and  a  reference  beam  (Ref)  as  de¬ 
picted  in  Figure  la.  in  this  case,  the  object  beam  is  shown  to  be  derived  from  an  auxilia¬ 
ry  optical  system.  During  the  raytrace  through  the  optical  system  under  consideration  (the 
"primary"  system),  a  reconstruction  ray  (Rec)  impinges  upon  the  hologram,  is  diffracted  by 
the  hologram  and  results  in  an  image  ray  (Img).  The  phase  and  direction  of  propagation  of 
the  image  ray  is  determined  by  the  phases  and  directions  of  the  reconstruction,  reference 
and  object  rays  at  the  reconstruction  ray  intercept.  A  simplified  form  of  the  grating 
equations  which  determine  the  image  ray  is  as  follows. 
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where  4  is  the  ray  phase,  )£  is  the  readout  wavelength,  is  the  recording  wavelength, 

1,  m  and  n  are  the  x,  y  and  z  direction  cosines,  respectively;  and  subscripts  I,  C,  0  and 
R  refer  to  the  image,  reconstruction,  object  and  reference  rays,  respectively.  The  sign 
choice  in  Equations  1-3  is  used  to  select  either  the  principal  diffracted  wavefront  (  +  )  or 
the  conjugate  wavefront  (-).  The  sign  choice  in  Equation  4  is  used  to  select  the  z-di- 
rection  of  propagation  of  the  wavefront. 


In  the  typical  raytrace  problem,  the  phase  and  direction  cosines  of  the  reconstruction 
ray  is  known,  either  as  input  to  the  system  or  as  the  result  of  a  raytrace  through  a  pre¬ 
ceding  element.  The  phases  and  direction  cosines  of  the  object  and  reference  rays,  how¬ 
ever,  must  be  determined  based  upon  the  reconstruction  ray  intercept.  For  the  case  of  a 
spherical  HOE  in  which  the  object  and  reference  wavefronts  are  restricted  to  being  either 
plane  or  spherical  wavef  :>nts,  the  task  is  a  simple  one.  The  phase  and  direction  cosines 
of  a  ray  passing  through  ny  point  illuminated  by  the  wavefront  are  easily  calculated  based 
upon  the  direction  of  the  plane  wavefront  or  the  location  of  the  spherical  point  source. 

The  case  of  an  aspheric  HOE,  however,  may  present  a  more  difficult  task  depending  upon  the 
way  in  which  the  arbitrary  recording  wavefront  is  defined. 

The  two  principle  ways  of  defining  an  arbitrary  wavefront  are  (1)  by  the  specification 
of  an  auxiliary  optical  system  that  is  used  to  generate  the  wavefront  (Figure  la)  or  (2) 
by  an  analytical  description  of  a  wavefront  defined  on  a  surface  (Figure  lb).  In  this  pa¬ 
per  we  discuss  both  ways  of  describing  arbitrary  recording  wavefronts,  with  emphasis  on 
the  analytical  description  on  a  surface. 

Arbitrary  recording  wavefronts  derived  from  auxiliary  optical  systems 

The  task  of  determining  the  phase  and  direction  cosines  of  a  recording  wavefront  derived 
from  an  auxiliary  optical  system  at  a  given  reconstruction  ray  intercept  is  indeed  a  dif¬ 
ficult  one.  There  are  two  basic  approaches  to  solving  this  problem.  The  first  approach 
is  to  trace  rays  through  the  auxiliary  system  in  an  iterative  fashion  until  the  ray  which 
passes  through  the  reconstruction  ray  intercept  is  found.  The  second  approach  is  to  trace 
a  grid  of  rays  through  the  auxiliary  system  to  the  HOE,  and  then  during  the  raytracing  of 
the  primary  system  to  perform  an  interpolation  on  the  grid  of  rays  to  obtain  the  phases 
and  direction  cosines  at  the  reconstruction  ray  intercepts.  A  two-dimensional  interpo¬ 
lation  is  required  which  preferably  takes  into  account  the  known  direction  cosine  samples 
as  well  as  the  phase  samples.  The  interpolation  is  further  complicated  by  the  fact  that  a 
regularly  spaced  grid  at  the  input  to  the  auxiliary  optical  system  will  be  distorted  into 
an  irregular  grid  at  the  HOE. 

Each  of  the  above  approaches  has  its  own  advantages  and  disadvantages.  The  iterative 
approach  requires  that  several  raytraces  through  the  auxiliary  system  be  performed  for  each 
ray  traced  through  the  primary  system.  In  general  the  number  of  iterative  raytraces 
through  the  auxiliary  system  for  each  primary  system  raytrace  will  be  small  since  the  arbi¬ 
trary  recording  wavefronts  of  interest  for  HOEs  are  well  behaved.  The  interpolation  ap¬ 
proach,  on  the  other  hand,  requires  that  a  relatively  large  number  of  rays  (typically  25 
to  several  hundred)  be  traced  through  the  auxiliary  system  once.  Thereafter,  any  number 


of  rays  may  be  traced  through  the  primary  system  and  only  interpolation  will  be  required 
to  determine  the  phase  and  direction  cosines  of  the  arbitrary  recording  wavefront  at  each 
reconstruction  ray  intercept.  If  the  auxiliary  system  contains  an  optimization  variable, 
however,  the  grid  of  rays  will  need  to  be  retraced  each  time  the  variable  changes  value. 

The  tradeoffs  between  the  two  approaches  are  complex  and  depend  in  large  part  on  the  rela¬ 
tive  complexities  of  the  primary  and  auxiliary  systems. 

The  interpolation  approach  has  been  implemented  within  the  Holographic  Optics  Analysis 
and  Design  (HOAD)2  program  at  ERIM  and  has  been  discussed  in  detail  elsewhere.2  The 
remainder  of  this  paper  will  concentrate  on  the  class  of  aspheric  HOEs  for  which  the  arbi¬ 
trary  recording  wavefronts  are  defined  analytically  at  the  recording  surface. 

Analytical  arbitrary  wavefronts  for  aspheric  HOE  design 

We  have  recently  added  to  the  HOAD  raytrace  program  at  ERIM  the  capability  of  analyzing 
aspheric  HOEs  recorded  with  analytically  defined  wavefronts.  The  analysis  is  currently 
limited  to  wavefronts  which  are  pure  phase  functions,  i.e.,  the  wavefront  amplitude  is  as¬ 
sumed  uniform  over  the  extent  of  the  wavefront.  The  phase  function  is  assumed  to  be  de¬ 
fined  on  ihe  surface  of  the  recording  medium.  The  surface  (or  substrate)  can  be  curved  as 
well  as  flat.  A  variety  of  general  analytical  phase  functions  have  been  supplied  for  the 
designer's  use  including  the  following. 

.  Sum  of  monomials  (power  series) 

9  9 

♦  <*.  y)  “I  I  C.  .  xjyJ 

i=0  j“0  lJ  <5) 

.  Sum  of  Legendre  polynomials  (orthogonal  polynomials) 

9  9 
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.  Spherical  wavefront  +  sum  of  monomials 
.  Spherical  wavefront  +  sum  of  polynomials 

The  Cij  in  Equations  (5)  and  (6)  represent  the  coefficient  values  of  the  respective  poly¬ 
nomials.  Up  to  100  coefficients  may  be  specified  for  each  phase  function  and  all  may  be 
used  as  optimization  parameters.  The  maximum  of  ninth  order  in  x  and  y  was  made  to  limit 
the  amount  of  computer  memory  required  to  store  the  coefficients  of  each  arbitrary  wave- 
front.  In  addition  to  using  the  preprogrammed  phase  functions  described  above,  the  de¬ 
signer  may  himself  define  an  explicit  phase  function  utilizing  up  to  100  optimization  para¬ 
meters.  New  phase  functions  are  continually  being  added  to  the  library  of  available 
functions  in  this  manner. 

There  are,  of  course,  limitations  upon  how  quickly  any  phase  function  may  vary.  Speci¬ 
fically,  in  order  to  avoid  evanescent  waves,  the  following  relation  must  be  satisfied, 

./S?  +  <  2&  <7) 

\»y )  -  * 

where  x  is  the  recording  wavelength.  Violation  of  this  constraint  results  in  a  ray  fail¬ 
ure  during  the  raytrace. 

Flat  substrate 

Raytracing  through  a  flat  aspheric  HOE  recorded  with  an  analytically  defined  arbitrary 
wavefront  is  a  straightforward  procedure.  The  direction  cosines  of  the  analytical  wave- 
front  at  the  reconstruction  ray  intercept,  (xo,  yo) •  ar*  readily  computed  from  the  par¬ 
tial  derivatives  of  the  phase  function  as  follows. 

♦  -  ♦(*,  y)L  „  (8) 
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(10) 


The  sign  of  the  z  direction  cosine  given  in  Equation  (11)  is  chosen  to  provide  the  desired 
direction  of  propagation.  Then,  using  the  grating  Equations  (1)  to  (4),  the  phase  and  di¬ 
rection  cosines  of  the  image  ray  are  computed. 

Curved  substrate 


The  use  of  an  analytically  defined  wavefront  to  record  an  aspheric  HOE  on  a  curved  sub¬ 
strate  is  somewhat  more  complex  than  the  flat  substrate  case.  One  method  of  specifying 
the  wavefront  would  be  to  define  a  phase  function  on  the  curved  surface,  in  which  case  the 
computation  of  the  direction  cosines  would  depend  upon  the  surface  function  in  addition  to 
the  phase  funtion  (Figure  3a) •  An  alternative  approach  would  be  to  define  the  phase  func¬ 
tion  on  a  plane  separated  from  the  curved  substrate.  This  latter  case  is  similar  to  that 
of  defining  a  wavefront  by  an  auxiliary  system  in  that,  in  order  to  find  the  ray  phase  and 
direction  cosines  at  a  given  intercept,  it  is  necessary  to  either  perform  an  iterative  ray- 
trace  (Figure  3b)  or  trace  a  grid  of  rays  and  use  iterpolation.  In  this  case  the  pro¬ 
cedure  is  simplified  by  the  fact  that  the  raytrace  is  between  two  surfaces  with  no  inter¬ 
vening  optics.  It  is  particularly  simple  if  the  plane  is  chosen  to  be  a  tangent  plane  of 
the  curved  surface. 


Figure  3.  A  curved  aspheric  HOE  recorded  with  an  analytically  defined  wavefront  can  have 
that  wavefront  (a)  defined  on  the  curved  surface  or  (b)  defined  on  a  plane  separated  from 

the  surface. 


He  have  chosen  to  implement  the  first  method  in  which  the  phase  function  is  defined  on 
the  curved  surface.  The  analyses  proceeds  as  follows.  Let  the  phase  of  a  wavefront 
throughout  a  certain  volume  be  4  (x,  y,  z) .  Assuming  that  the  wavefront  forms  a  normal 
congruence  (i.e.,  rays  do  not  cross  one  another  and  the  phase  and  direction  cosines  are 
uniquely  defined  at  every  point),  the  wavevector  is  given  by 

C  -  kxx  +  kyy  +  kzz  -  vj(x.  y,  z)  (12) 

where  x,  y,  and  z  are  the  cartesian  unit  vectors.  The  direction  cosines  are  given  by 
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whece  Xq  is  the  wavelength  of  light  in  a  vacuum  and  nj  is  the  index  of  refraction.  The 
phase  evaluated  at  a  surface 

z  -  z (x,  y)  (17) 

is  given  by 

♦(*.  y)  “  ♦!*.  y.  z(x.  y) 1  (18) 

which  we  refer  to  as  the  surface  phase  function.  Having  defined  the  surface  phase  function 
♦  (x,  y)  by  an  analytical  expression,  the  problem  is  to  recover  the  direction  cosines  which 
are  proportional  to  the  partial  derivatives  of  * (x,  y,  x) ,  which  are  different  fro*  the 
partial  derivatives  of  +(x,  y) .  By  the  chain  rule  of  partial  derivatives,  we  have  the  two 
equations 

It  -  *  _  it  +  it  IS.  (19) 
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and,  in  addition,  taking  the  magnitude  squared  of  Equation  (12)  we  have  the  condition 
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Solving  the  three  simultaneous  Equations  (19),  (20)  and  (21)  for  the  three  unknowns,  we 
find  that 
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where  zx  »  az(x,  y)/ax  and  Zy  >  az(x,  y)/ay.  The  computation  of  the  phase  and  direc¬ 
tion  cosines  from  the  surface  phase  function  using  Equations  22  to  24  is  very  straightfor¬ 
ward  in  contrast  with  the  method  of  defining  the  phase  on  a  plane  separated  from  the  curved 
surface.  Unlike  that  method,  however,  the  definition  of  the  phase  on  the  curved  surface 
causes  a  strong  coupling  between  the  wavefront  +  (x,  y,  z)  and  the  shape  of  the  surface 
z(x,  y). 

Recording  techniques  for  aspheric  HOBs 

Once  a  recording  wavefront  is  established  by  the  designer,  it  becomes  necessary  to  pro¬ 
duce  that  wavefront  in  the  laboratory.  For  the  case  of  a  recording  wavefront  defined  by 
an  auxiliary  system  as  depicted  in  Figure  la,  the  procedure  is  straightforward:  an  optical 
system  corresponding  to  the  auxiliary  system  must  be  assembled.  For  the  case  of  an 
analytically  defined  arbitrary  wavefront,  on  the  other  hand,  the  design  itself  does  not 
suggest  a  method  of  arriving  at  the  desired  wavefront. 

One  method  would  be  to  record  the  desired  aspheric  HOE  as  a  CGH.  Several  inherent  limi¬ 
tations  of  CCHs,  however,  severely  limit  the  direct  use  of  a  CGH  as  a  HOE.  First  and  fore¬ 
most,  optical  recording  devices  used  to  generate  CGHs  are  limited  in  spatial  resolution 
and  space-bandwidth  product,  restricting  the  angles  of  diffraction  and  numerical  apertures 
of  the  CGH.  in  addition,  unwanted  orders  of  diffraction  and  other  spurious  terms  are  usu¬ 
ally  present  in  a  CGH.  These  further  restrict  the  usable  diffraction  angles  and 
f ield-of-view  If  interference  with  these  undesired  terms  is  to  be  avoided.  Another  limi- 


tation  is  the  low  value  of  the  maximum  diffraction  efficiencies  for  most  types  of  CGHs. 

The  maximum  diffraction  efficiency  for  a  thin  amplitude  hologram  is  6.25*,  for  a  binary 
amplitude  hologram  10.13%,  and  for  a  thin  phase  hologram  33.9*. 5  Types  of  CGHs  that 
have  diffraction  efficiencies  approaching  100%,  such  as  the  kinoform®  and  the  ROACH,7 
are  more  difficult  to  generate  accurately. 

All  of  the  above  limitations  can  be  circumvented,  however,  through  the  use  of  a  HOE  re¬ 
corded  with  a  CGH  in  one  of  the  recording  beams,  instead  of  using  the  CGH  itself  as  the 
optical  element.  By  using  a  volume  phase  material,  such  as  dichromated  gelatin®  for  the 
HOE,  diffraction  efficiencies  approaching  100*  are  achieved  and  spurious  orders  of  diffrac 
tion  are  minimized.  By  combining  appropriate  optics  with  the  CGH  in  the  recording  beam, 
the  resolution  and  the  space-bandwidth  product  required  of  the  CGH  can  be  greatly 
reduced.’  In  addition,  spatial  filtering  can  be  performed  in  order  to  remove  spurious 
terms  inherent  in  the  CGH. 

The  first  step  is  to  design  a  recording  system  that  reduces  the  space-bandwidth  product 
and  resolution  requirements  of  the  CGH.  It  is  assumed  that  the  recording  wavefront  of  the 
HOE  has  already  been  specified  as  an  analytical  arbitrary  wavefront.  Then  the  raytrace 
program  is  used  to  back  propagate  the  desired  recording  wavefront  through  a  recording  op¬ 
tical  system  to  a  CGH  plane  (Figure  4a) .  The  recording  optical  system  will  ordinarily  be 
designed  to  remove  tilt,  focus,  and  other  low-order  phase  terms  which  tend  to  be  large  in 
magnitude.  A  secondary  purpose  of  the  recording  optical  system  is  to  provide  a  frequency 
plane  in  which  a  spatial  filter  may  be  used  to  remove  the  undesired  diffracted  orders  of 
the  CGH.  Once  the  recording  optical  system  is  designed  to  produce  a  wavefront  with  accept 
ably  low  space-bandwidth  product  at  the  CGH  plane,  a  grid  of  rays  is  back  propagated  from 
the  hologram  to  the  CGH  to  provide  samples  of  the  phase  fucntion  that  is  to  be  recorded  as 
a  CGH.  A  frequency  offset  must  be  added  during  the  CGH  recording  process  to  insure  that 
the  zero-order  and  second-order  diffracted  terms  do  not  overlap  the  desired  first  order 
diffracted  term  in  the  frequency  plane.  The  amplitude  transmittance  of  the  CGH  then  has 
the  form 

t  “  b  +  m  cos  [ax  +  4>(x,  y)  ]  (25) 

where  b  =0.5  is  a  bias,  m  <_  0.5  is  the  modulation,  and  ux  is  the  carrier  frequency  offset. 

For  this  Burch-type10  CGH  (simple  carrier  frequency) ,  an  offset  slightly  larger  than  one 
half  of  the  double-sided  bandwidth  of  the  wavefront  is  required.  This  is  less  than  that 
generally  required  by  an  optically  generated  hologram  because  the  Burch  type  CGH  does  not 
record  the  object  autocorrelation  term  which  has  twice  the  bandwidth  of  the  object  wave- 
front.  However,  in  order  to  avoid  spurious  terms  that  would  arise  in  the  event  of  a  non¬ 
linear  amplitude  transmittance,  it  may  be  necessary  to  use  a  carrier  frequency  that  is  1.5 
times  the  double-sided  bandwidth  of  the  wavefront.’ 


CGH  Plane 
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Figure  4.  COHOE  recording  process,  (a)  After  the  desired  arbitrary  recording  wavefront 
Is  defined  It  Is  backwards  raytraced  through  an  auxiliary  recording  system  to  a  CGH  defi¬ 
nition  plane,  (b)  The  COHOE  is  recorded  on  a  high  efficiency  medium  using  the  CGH  and 

optical  system  as  defined  in  (a) . 


Finally,  the  CGH  is  fabricated  and  assembled  with  the  recording  system  (Figure  4b)  to 
provide  the  desired  recording  wavefront  at  the  HOE,  The  HOE  which  is  recorded  in  this  man¬ 
ner  we  call  a  computer-originated  HOE  (COHOE).  In  the  following  sections  we  discuss  an 
example  of  the  design  of  an  aspheric  HOE  and  its  implementation  as  a  COHOE, 

An  aspheric  HOE  design  example 

To  demonstrate  the  use  of  an  analytical  arbitrary  recording  wavefront  we  designed  a 
Fourier  transform  HOE  which  could  be  used  in  a  coherent  optical  processor  (Figure  5).  A 
transparency  at  the  input  plane  is  illuminated  by  a  coherent  plane  wavefront.  The  input 
transparency  produces  an  angular  spectrum  of  plane  wavefronts  (one  for  each  spatial  fre¬ 
quency  component  of  the  input)  which  propagate  to  the  Fourier  transform  HOE.  The  HOB 
causes  them  to  be  focused  to  points  in  the  output  plane.  The  higher  spatial  frequencies 
in  the  input  transparency  diffract  the  illuminating  wavefront  at  proportionately  higher 
angles  and  come  to  focus  farther  from  the  center  of  the  output  plane. 


Figure  5.  Fourier  transform  HOE  readout  geometry.  Three  bundles  of  five  parallel  input 
rays  each  are  shown  propagating  from  the  input  plane  to  the  Fourier  transform  plane.  Each 
of  the  three  ray  bundles  represents  a  different  plane-wave  spatial  frequency  component  which 
comes  to  focus  at  a  point  In  the  transform  plane. 

In  a  previous  design  effort  to  produce  such  a  Fourier  transform  HOE11  using  conven¬ 
tional  spherical  wavefronts,  it  was  determined  that  optimum  performance  over  a  range  of 
input  spatial  frequencies  was  achieved  with  a  recording  geometry  as  shown  in  Figure  6. 

The  point  source  for  the  object  recording  wavefront  should  be  on  an  axis  normal  to  the  HOE 
at  a  point  corresponding  to  the  center  of  the  output  plane. 


Figure  6.  The  optimized  recording  geometry  for  a  conventional  Fourier  transform  HOE  using 

plane  and  spherical  wavefronts. 


Using  this  previous  design  as  a  starting  point,  we  further  optimized  the  HOB,  allowing 
the  tilted  plane  reference  wavefront  to  be  perturbed  by  the  following  polynomial  phase 
function. 
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All  ten  of  the  Ctj  coefficients  were  allowed  to  vary  during  a  damped  least-squares 
optimization.  Twenty-one  rays  forming  a  pair  of  orthogonal  fans  on  a  25  mm  diameter  input 
aperture  (Figure  7a)  were  raytraced  to  the  Fourier  transform  plane  for  each  system  solution 
in  the  optimization.  The  focal  length  of  the  HOE  was  designed  to  be  0.5  meters  and  the  re¬ 
cording  and  readout  wavelengths  were  both  514.5  nra.  The  merit  function  consisted  of  the 
sum  of  squares  of  the  RMS  spot  size  at  the  Fourier  transform  plane  for  ten  illumination 
angles:  a  *  -2.4°,  -1.2°,  0°,  1.2°,  2.4°  (in  the  x-z  plane)  and  8  «  -2.4°, 

-1.2°,  0°,  1.2°,  2.4°  (in  the  y-z  plane).  The  resultant  optimized  coefficients 
are  given  in  Table  1.  The  coefficients  of  Table  1  are  normalized  to  have  units  of  wave¬ 
lengths  (X)  and  the  x  and  y  coordinates  are  scaled  such  that  -1  i  x,  y  <  1  over  the  holo¬ 
gram  recording  area  (i.e.,  x  and  y  are  unitless).  A  perspective  plot  of  the  optimized 
phase  correction  is  shown  in  Figure  8.  Aberration  plots  which  compare  the  primary  aber¬ 
rations  of  the  starting  design  (conventional  spherical  HOE)  and  the  optimized  aspheric  HOE 
as  a  function  of  field  angle  a (x-z  plane)  are  depicted  in  Figure  9. 


(c)  Spherical 


(d)  Total  RMS 


Figure  9.  Aberration  plota  comparing  the  performance  of  a  conventional  Fourier  transform 
HOE  (x's)  and  an  aspheric  Fourier  transform  HOE  (+'s). 


Recording  and  evaluation  of  an  aspheric  HOE 


The  aspheric  Fourier  transform  HOE  design  described  in  the  last  section  was  recorded  as 
a  COHOE  using  the  optical  system  diagrammed  in  Figure  10.  The  diffracted  wavefront  from 
the  CGH  is  reimaged  at  the  COHOE  recording  plane  by  a  one-to-one  telescope.  This  telescope 
not  only  performs  imaging  from  the  CGH  plane  to  the  COHOE  plane,  but  it  also  preserves  the 
desired  phase  relationships  (i.e.,  it  does  not  introduce  an  extra  spherical  phase  term  as 
would  imaging  with  a  single  lens).  For  this  case  it  was  assumed  that  the  imaging  system 
adds  no  extra  phase  terms,  so  in  the  design  it  was  not  necessary  to  simulate  the  effects 
of  the  imaging  system.  A  spatial  filter  mask  is  positioned  at  the  frequency  plane  such 
that  only  the  desired  first-order  diffracted  wavefront  of  the  CGH  is  passed  to  the  record¬ 
ing  plane.  The  tilt  of  the  COHOE  is  such  that  the  desired  20°  offset  angle  is  obtained 
at  the  recording  plane.  A  slight  additional  tilt  is  added  to  allow  for  the  angle  that  the 
first  order  diffracted  term  makes  with  the  optical  axis  as  it  exits  from  the  telescope. 

The  CGH  is  similarly  tilted  so  that  it  and  the  COHOE  are  in  conjugate  image  planes.  An 
objective  and  pinhole  assembly  provides  the  required  point  source  object  beam. 


Figure  10.  The  COHOE  recording  geometry  used  to  produce  a  Fourier  transform  HOE  with  an 
aspheric  reference  wavefront  as  defined  by  Table  1. 

The  optimized  arbitrary  recording  wavefront  described  by  Equation  (25)  and  Table  1  has 
maximum  bandwidths  of  0.747  cyc/mm  in  the  x  dimension  and  0.758  cyc/mm  in  the  y  dimension 
when  evaluated  over  a  90  mm  diameter  aperture.  This  wavefront  was  recorded  as  a  CGH  on  an 
Optronics  Model  1600  film  recorder  after  a  carrier  frequency  of  3  cyc/mm  and  a  bias  were 
added.  The  film  recorder  was  operated  with  a  50  urn  square  recording  spot  on  a  50  urn  sample 
spacing.  The  Optronics  Model  1600  film  recorder  is  also  capable  of  recording  a  25  ym  aper¬ 
ture  on  a  25  m  sample  spacing  with  reduced  speed.  The  Optronics  film  recorder  has  a  ro¬ 
tating  drum  with  a  translating  LED  and  uses  Kodak  Linagraph  Shellburst  film  #2474.  The 
CGH  was  contact  copied  onto  a  Kodak  #649F  microflat  plate  for  insertion  into  the  COHOE  re¬ 
cording  system.  The  microflat  substrate  minimized  undesired  phase  errors  due  to  variations 
in  substrate  thickness.  Similarly,  the  COHOE  was  recorded  on  a  Kodak  #131-01  microflat 
plate  to  insure  that  phase  errors  were  not  introduced  during  readout.  The  recording  of 
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Figure  11.  Interferograms  of  the  aspheric  wavefront  defined  by  Table  1  (a)  as  predicted 
by  a  computer  reytrace  end  (b)  as  optically  recorded. 


the  COHOE  utilized  514.5  iui  light  from  an  Argon-ion  laser.  A  computer-predicted  interfero- 
graa  of  the  desired  wavefront  correction  was  plotted  (Figure  11a)  and  compared  with  an  op¬ 
tically  derived  interf erogram  of  the  wavefront  as  recorded  by  the  CO HOE  (Figure  lib).  The 
agreement  was  excellent,  indicating  that  no  significant  phase  errors  were  inadvertently 
introduced  by  the  recording  system. 

The  performance  of  the  COHOE  was  evaluated  using  a  rotatable  mirror  positioned  at  the 
input  transparency  plane  (Figure  12).  in  this  manner,  a  single  stationary  plane  wavefront 
was  made  to  simulate  the  wavefront  that  would  be  produced  by  any  single  spatial  frequency 
at  the  input  transparency  plane.  A  microscope  attached  to  a  precision  translation  device 
was  positioned  such  that  it  was  focused  on  the  Fourier  transform  plane  of  the  COHOE. 


Figure  12.  Optical  arrangement  for  evaluating  the  performance  of  the  aspheric  Fourier 
transform  COHOE  using  a  rotatable  mirror  (M2)  at  the  input  transparency  plane  to  simulate 
any  single  plane-wave  spatial  frequency  cosiponent. 

The  rotatable  mirror  was  positioned  sequentially  at  the  design  field  angles  and  the  re¬ 
sulting  spot  sizes  in  the  Fourier  transform  plane  were  recorded  on  film.  A  comparison  of 
the  computer  predicted  spot  sizes  of  Figure  13  with  the  corresponding  measured  spot  sizes 
of  Figure  14  indicates  excellent  agreement.  The  computer  generated  spot  diagrams  were  gen¬ 
erated  by  raytracing  a  hexapolar  distribution  of  rays  (Figure  7b)  at  each  of  the  specified 
field  angles.  These  spot  measurements  were  made  with  an  input  aperture  size  of  35.6  mm 
which  is  somewhat  larger  than  the  design  size  of  25.4  mm.  The  impulse  response  of  the 
COHOE  recording  system  has  been  included  in  Figure  14  for  comparison.  The  predicted  im¬ 
provement  in  performance  of  the  COHOE  as  compared  with  the  conventional  (spherical)  HOE, 
particularly  at  the  large  field  angles,  is  verified. 
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Figure  13.  Computer  predicted  spot  sizes  of  the  conventional  and  computer  originated 

Fourier  transform  HOEs. 
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Figure  14.  Optically  recorded  spot  sizes  of  the  conventional  and  computer  originated 
Fourier  transform  HOEs  using  the  setup  shown  in  Figure  12. 


Conclusions 


In  summary,  the  work  reported  In  this  paper  has  demonstrated  the  feasibility  of  analyz¬ 
ing  and  Implementing  aspheric  holographic  optical  elements  using  analytic  descriptions  of 
the  recording  wavefronts  during  the  design  phase  and  using  a  CGH  in  the  recording  beam  dur¬ 
ing  fabrication.  The  inherent  properties  of  computer-generated  holograms  which  limit  their 
direct  use  as  holographic  optical  elements  are  avoided  by  using  the  COHOE  recording  tech¬ 
nique.  An  aspheric  Fourier  transform  HOE  was  designed  which  demonstrated  a  significant  im¬ 
provement  in  performance  when  compared  to  a  conventional  HOE  recorded  with  spherical  wave- 
fronts.  It  is  our  opinion  that  aspheric  HOEs  will  prove  to  be  generally  valuable  in  the 
design  of  future  high-performance  holographic  optical  systems.  In  particular  we  expect 
that  aspheric  HOEs  will  provide  better  performance  with  fewer  elements  and  will  lessen  the 
burden  on  the  refractive  optics  in  hybrid  systems. 
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Appendix  B 


TITLE:  Checkerboard  Real -Imaginary  Phase  Code 

AUTHOR:  J.R.  Fienup 

ABSTRACT:  This  phase  code  (diffuser)  for  holography  causes  the  spectrum  to  be 
leveled  well  (reduced  dynamic  range)  and  results  In  a  minimum  of 
speckle  in  the  Image. 

SUMMARY: 

In  digital  holography,  phase  codes  (which  are  analogous  to  diffusers 
In  Interferometric  holography)  are  applied  to  the  object  In  order  to 
reduce  the  dynamic  range  of  Its  Fourier  transform,  which  is  encoded  In 
the  computer-generated  hologram.  This  makes  it  easier  to  faithfully 
record  the  hologram  and  Increases  Its  diffraction  efficiency.  An  unde¬ 
sirable  artifact  arising  from  most  phase  codes  (particularly  the  random 
phase  code)  Is  speckle  In  the  reconstructed  image.  A  smoothed,  slowly 
varying  phase  code  results  In  less  speckle  in  the  Image  but  a  greater 
dynamic  range  in  Its  Fourier  transform.  The  Checkerboard  Real -Imaginary 
Phase  (CRIP)  code,  when  combined  with  an  appropriate  weighting  at  the 
hologram  aperture,  results  In  both  a  greatly  reduced  dynamic  range  and 
greatly  reduced  speckle.  The  CRIP  code  consists  of  a  checkerboard  pat¬ 
tern  of  +1  (pure  real)  and  ±1  (pure  Imaginary),  where  the  signs  are  chosen 
quasi-randomly.  A  physical  Interpretation  of  the  origin  of  the  speckle 
will  be  described,  showing  why  the  CRIP  code  results  In  greatly  reduced 
speckle.  Experimental  results  will  be  shown  comparing  the  CRIP 
code  with  other  phase  codes. 


Published  In: 

J.  Opt.  Soc.  Am.  68,  p.  1444  (October  1978). 


